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Nematicity and magnetism are two key features in Fe-based superconductors, and their interplay 
is one of the most important unsolved problems. In FeSe, the magnetic order is absent below the 
structural transition temperature Tatr = 90K, in stark contrast that the magnetism emerges slightly 
below Tstr in other families. To understand such amazing material dependence, we investigate the 
spin-fluctuation-mediated orbital order (uxz ^ nyz) by focusing on the orbital-spin interplay driven 
by the strong-coupling effect, called the vertex correction. This orbital-spin interplay is very strong 
in FeSe because of the small ratio between the Hund’s and Coulomb interactions (J/tJ) and large 
dxz, dyz-orhitaXs weight at the Fermi level. For this reason, in the FeSe model, the orbital order 
is established irrespective that the spin fluctuations are very weak, so the magnetism is absent 
below Tstr. In contrast, in the LaFeAsO model, the magnetic order appears just below Tstr both 
experimentally and theoretically. Thus, the orbital-spin interplay due to the vertex correction is the 
key ingredient in understanding the rich phase diagram with nematicity and magnetism in Fe-based 
superconductors in a unified way. 

PACS numbers: 74.70.Xa, 75.25.Dk, 74.20.Pq 


I. INTRODUCTION 

In Fe-based superconductors, the origin of the elec¬ 
tronic nematic state and its relation to the magnetism 
have been a central unsolved problem. Recently, the non¬ 
magnetic nematic state in FeSe has attracted increasing 
attention as a key to solve the origin of the nematic¬ 
ity. FeSe undergoes a structural and superconducting 
transitions at Tstr = 90K and Tc = 9K, respectively, 
whereas the magnetic transition is absent down to 0 K 
[1]. The strength of the low-energy antiferro-magnetic 
(AFM) fluctuations is very weak above Tgtr, while it 
starts to increase below Tgtr [2-7]. In stark contrast, the 
magnetic transition occurs at Tmag slightly below Tstr 
in other undoped Fe-based superconductors. Since the 
relation Tstr > Tmag is unable to be explained by the 
random-phase-approximation (RPA), we should develop 
the microscopic theory beyond the mean-field-level ap¬ 
proximations. 

Up to now, two promising triggers for the structure 
transition have been discussed intensively: In the spin- 
nematic scenario [8-12], the trigger is the spin-nematic 
order. This spin-fluctuation induced spin-quadrupole or¬ 
der could emerge above Tmag in highly magnetically frus¬ 
trated systems. In the orbital order scenario [13-16], 
the trigger is the ferro-orbital (FO) order nxz ^ nyz. 
Above Tstr, the strong orbital or spin-nematic fluctua¬ 
tions are observed by the measurements of shear mod¬ 
ulus Cee [2, 17, 18], Raman spectroscopy [19-22], and 
in-plane resistivity anisotropy [23, 24]. The nematic or¬ 
bital fluctuations originate from the strong orbital-spin 
mode-coupling due to the strong-coupling effect, which 
is described by the Aslamazov-Larkin vertex correction 
(AL-VC). The electronic nematic state studied in single¬ 
orbital models [25] is more easily realized in multiorbital 
systems thanks to the AL-VC mechanism [16]. 


Except for the presence or absence of magnetism below 
Tstr, FeSe and other Fe-based superconductors show com¬ 
mon electronic properties. Below Tgtr, in both FeSe and 
BaFe 2 As 2 , large orbital polarization AE = Eyz — E^z ^ 
50 meV [26-33] is observed. Such large AE originates 
from the electron-electron correlation since the lattice 
distortion {a — b)/{a-\- b) is just 0.2 ~ 0.3%, as we dis¬ 
cuss based on band calculation in Appendix A. Above 
Tstr, the electronic nematic susceptibility is enhanced in 
both BaFe 2 As 2 [17, 19, 23] and FeSe [2, 18], following the 
similar Curie-Weiss behavior. These facts indicate that 
the common microscopic mechanism drives the nematic 
order and fluctuations in all Fe-based superconductors, 
in spite of the presence or absence of the magnetism. 

The realistic multiorbital Hubbard models for Fe-based 
superconductors, which are indispensable for the present 
study, were derived by using the first-principles method 
in Ref. [34]. To understand the absence of the magnetism 
below Tstr in FeSe, one significant hint is the smallness of 
the ratio between the Hund’s and Coulomb interactions, 
J/U, since the Hund’s coupling enlarges (suppresses) the 
intra-site magnetic (orbital) polarization, which is veri¬ 
fied by the functional renormalization-group (fRG) the¬ 
ory [35, 36]. Another significant hint is the absence of the 
(ia;j,-orbital hole-pocket in FeSe, which is favorable for the 
orbital-spin interplay on the {d^z, dyz)-orhitdls due to the 
AL-VC mechanism. 

The goal of this paper is to explain the amazing variety 
of the electronic nematic states in Fe-based superconduc¬ 
tors, especially the non-magnetic nematic state in FeSe, 
on the same footing microscopically. For this purpose, 
we study the spin-fluctuation-mediated orbital order by 
applying the self-consistent vertex-correction (SC-VC) 
method [16] to the first-principles models. In FeSe, the 
orbital-spin interplay is significant because of the small¬ 
ness of J/U and the absence of dxy-\io\e pocket. For 


2 


this reason, the orbital order is realized even when the 
spin fluctuations are substantially weak. The rich vari¬ 
ety of the phase diagrams in Fe-based superconductors, 
such as the presence or absence of the magnetic order 
in the nematic phase, are well understood by analyzing 
the vertex correction seriously. The SC-VC theory had 
been successfully applied to explain the phase diagram 
in LaFeAs(0,H) [37], nematic CDW in cuprates [38, 39], 
and triplet superconductivity in Sr 2 Ru 04 [36]. 

We comment that the localized spin models have been 
successfully applied to the nematic order, stripe magnetic 
order, and so on [40]. On the other hand, weak-coupling 
theories have also been applied to Fe-based superconduc¬ 
tors satisfactorily [41]. In the present study, we study the 
mechanisms of the nematicity and magnetism in various 
Fe-based superconductors in terms of the itinerant pic¬ 
ture, by taking the strong-coupling effect due to the AL- 
VC into account. The significant role of the AL-VC on 
the orbital fluctuations has been confirmed by the fRG 
theory [35, 36]. The AL-VC is important to reproduce 
the Kugel-Khomskii-type orbital-spin interaction [16]. 


II. MODEL HAMILTONIAN AND SC-VC 
THEORY 


shallow, the superconductivity in FeSe could be close to 
a BCS-BEC crossover [45]. 

In Eq. (1), is the first-principles screened Coulomb 
potential for d-orbitals given by the “constrained-RPA 
method” [34] given as 




where Um,i and Jm,i are orbital-dependent Coulomb 
and Bund’s interactions for d-electrons, respectively [34]. 
The averaged intra-orbital Coulomb interaction U = 
i Ya=i UiJ Bund’s interactions J = ^ J2i>m 
are (C/, J)=(7.21, 0.681) for FeSe, and ([/, J)=(4.23, 
0.568) for LaFeAsO in unit eV [34]. Thus, the ratio 
J/U = 0.0945 in FeSe is much smaller than the ratio 
J/U = 0.134 in LaFeAsO. Such strong material depen¬ 
dence of {U, J) is understood as follows: Ui^m is strongly 
screened by the screening bands (excluding the 8 bands in 
iL^) whereas the screening of Ji^m is much weak, and the 
number of the screening bands is small in FeSe [46]. The 
factor r(< 1) in Eq. (1) is introduced to adjust the spin 
fluctuation strength. The ratio Ji,m/Ui^rn is unchanged 
by introducing the factor r [47, 48]. 


In the present study, we study the realistic d-p Bub- 
bard models 


Huir) — 77 m + ( 1 ) 

for M=LaFeAsO and FeSe by applying the SC-VC 
method [16]. In Eq. (1), 

~ ( 2 ) 

k,lm,(7 

is the 8-orbital d-p tight-binding (TB) model in fc-space, 
which is obtained by using the WIEN2k and WAN- 
NIER90 softwares; see Appendix A for detailed expla¬ 
nation. a is the spin index, and l,m are the orbital in¬ 
dices: Bereafter, we denote the five d-orbitals as 
dxz, dyz, dxy, dx^-y 2 as 1,2,3,4, 5, and three p-orbitals 
as 6 ~ 8. The bandstructure and Fermi surfaces (FSs) in 
the LaFeAsO model are shown in Figs. 1 (a) and (b), re¬ 
spectively. Similar FSs with three hole-like FSs (h-FSs) 
and two electron-like FSs (e-FSs) exist in many Fe-based 
superconductors. In FeSe, however, h-FS3 composed of 
d 2 ;y-orbital is absent, and the size of each FS is very 
small as clarified by the ARPES [27, 32, 33] and dBvA 
[42, 43] studies. To reproduce experimental bandstruc¬ 
ture of EeSe, we introduce the additional intra-orbital 
hopping parameters into TlpeSei order to shift the dxy- 
orbital band [d^^^/y^-orbital band] at (F, M, X) points by 
(0, —0.25, -1-0.24) [(—0.24, 0, +0.12)] in unit eV; see Ap¬ 
pendix A. These energy shifts might be induced by the 
self-energy [44]. The constructed FSs in the FeSe model 
is shown in Fig. 1 (c). Since each Fermi pocket is very 
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FIG. 1: (color online) (a) Bandstructures of the eight-orbital 
TB models for LaFeAsO and FeSe. (b) FSs for the LaFeAsO 
TB model, (c) FSs for the FeSe TB model. The colors corre¬ 
spond to 2 (green), 3 (red), and 4 (blue), respectively. 


The 8x8 Creen function in the orbital basis is given 
as 

G{k) = (z“^fe„ + M 


( 4 ) 
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where k = {k,en = (2n + l)TrT), h^{k) is the kinetic 
term in Eq. (2), and z~^ = 1 — cJS/9e|e=o represents the 
mass-enhancement due to the self-energy at the Fermi 
level. Here, we introduce the constant mass-enhancement 
factor for d-orbital l/zi{> 1). Then, Eq. (4) gives the 
coherent part of the Green function, which mainly deter¬ 
mines the low-energy electronic properties. In the present 
study, r and zi are the fitting parameters. In FeSe, 
the orbital order is obtained in the real first-principles 
Hamiltonian (r Ri 1) by taking the experimental mass- 
enhancement factors Ri 4 into account, as shown 
later. 

The d-orbital charge (spin) susceptibilities (per spin) 
is given in the following 5^ x 5^ matrix form: 

X‘=(")(q) = l>'=(")(q)(l - H«)d^('<)(q))-i (5) 

where <l>'^('’^(q) = x°(g) -|- is the irreducible sus¬ 

ceptibility for the charge (spin) channel. In the SC-VC 
theory, we employ the AL-VC as and perform 

the self-consistent calculation with respect to the AL-VC 
and susceptibilities. Using the Green function in Eq. (4), 
the bare susceptibility is 

XM';m.m'(9) = -TGl^m{k + q)Gm',1'{k), (6) 

k 

where q = (g,w; = 21ttT). Also, the AL-VC for the 
charge susceptibility is given as 

^l,l’-,a,b-ej{q]p) 

p ar^h 

+ QWeJ-g,hi~P) + ^a,b-,c,d{P + 9We,f-,g,hi~P)} 

^ ■^m,m';c,d;g,h{^^ P') ^ (^) 

where p = (p, Wm), and U®.'=(g) = -|- f®’'=x®’‘^(g)f®’'=. 

The three-point vertex A{q]p), which gives the coupling 
between two-magnon and one-orbiton, is given as 

Ald'-^a,b;e,f {q^P) 

= -Tj2GiAk' + q)Gfp{k')GbAk' - P), ( 8 ) 

k' 

and AAm'-c,d-,g,hi^''P) = Ac^h-m,g-,d,m'Ap) + 

Ag^d-m,c-h,m'iT,-P - q)- We stress that the strong 
temperature dependence of the three-point vertex is 
significant for realizing the orbital order. We include all 
I/^-terms without the double counting in order to obtain 
quantitatively reliable results. Equation (7) means that 
the charge AL-VC becomes large in the presence of 
strong spin fluctuations. More detailed explanations are 
presented in the textbook [49] 

In the present study, we neglected the spin-channel 
VCs since it is expected to be unimportant as discussed 
in Ref. [50]. In Appendix B, we verify validity of 
this simplification in the present model by performing a 
time-consuming self-consistent calculation with respect 


to both charge- and spin-channel Maki-Thompson (MT) 
and AL-VCs. 

Hereafter, we mainly discuss the total spin susceptibil¬ 
ity, X"{q) = and the orbital suscepti¬ 

bilities for 0 ^ 2 _y 2 = n^z-ny^, Xl-^-y-^{q) = X 2 , 2 ; 2 ; 2 ( 9 ) + 
X3.3;3.3('7) - 2xi.2;3.3(9)- The divergence of x^ 2 _j^ 2 (q) 
at q = 0 gives rise to the FO order Uxz nyz- The 
charge (spin) Stoner factor ac{s) is given by the max¬ 
imum eigenvalue of (q) in Eq. (5), and the 

charge (spin) susceptibility is enlarged in proportion to 
the charge (spin) Stoner enhancement factor Sc(s) = 
(1 - ac(S))“^- 

As explained in Ref. [50], the development of 
Xx^-y^iO) is mainly induced by the diagonal elements 

of with respect to / = 2,3. If we drop the off-diagonal 
elements of approximately, Xx^-y^iO) is given as 

X= 2 _, 2 (0) « 2$7(I - (I - 5J/C/)C/4>'=), (9) 

where U = C/ 2,2 = C/ 3 , 3 , J = A, 3 , = Xu;i,j(0) + 

^ii i z(^) (^ = 2 or 3). Thus, the charge Stoner factor is 
ac = (I —5J/C/)C/$'^ PS {1 — 5J/U){1 + UX’^), considering 
the relation x°(< 7 ) « I/C/. Within the RPA ($= = $* = 
X*^), only the spin fluctuations develop since the relation 
as > ac is satisfied for J > 0. However, the opposite 
relation ac > as is realized if the relation ^ x° is 
satisfied due to the charge-channel AL-VC [37]. 

III. NUMERICAL RESULTS FOR LAFEASO 
AND FESE 

First, we analyze the LaFeAsO model based on the 
SC-VC theory. For z = I for each I, the obtained x*(9) 
and Xx^-y 2 {Q) are shown in Fig. 2 (a) and (b), respec¬ 
tively, for r = 0.41 [U = 1.74 eV) at T = 50 meV. Here, 
the number of fe-meshes is 32 x 32, and the number of 
Matsubara frequencies is 256. Thus, both AFM and FO 
susceptibilities develop divergently, and the realized en¬ 
hancement factors are S's ~ 40 and Sc ~ 50. The r- 
dependences of the enhancement factors at T = 50 meV 
are shown in the inset of Fig. 2 (c): Both Ss and Sc 
increase with r, and they are equivalent at r* = 0.41. 
The lower the temperature is, the smaller r* is, whereas 
the value of Ss = Sc at r* is approximately independent 
of T. Similar result is obtained in BaFe 2 As 2 model as 
shown in Appendix C. 

In addition, other antiferro-orbital susceptibilities 

XlziQ) = 2 [x§. 4 ; 3 , 4 (g) + X 3 , 4 ; 4 . 3 (g)] and xlzio) = 

2 [X 2 , 4 ; 2 , 4 ( 9 ) + X 2 , 4 ; 4 , 2 ( 9 )] develop secondary as reported 
in previous studies [16, 37, 49]. The obtained results are 
essentially similar to the results obtained in the five d- 
orbital Hubbard model for LaFeAsO explained in Ref. 
[16]. 

Figure 2 (c) shows the temperature dependences of 
the Stoner enhancement factors at r = 0.41. Both 
Sc and Ss approximately follow the Curie-Weiss be¬ 
haviors with the charge and spin Weiss temperatures 
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6c = 48 meV and 9$ = 40 meV, respectively. The ob¬ 
tained relation Xx^-y 2 ( 0 ) oc {T — 9c)~^ is consistent with 
the Curie-Weiss behavior of the nematic susceptibilities 
in BaFe 2 As 2 , derived from Cqq [17, 18], Raman spec¬ 
troscopy [19, 20], and in-plane resistivity anisotropy [23]. 
Since 0c 6s, one could interpret that the orbital order 
in LaFeAsO is driven by the spin fluctuations. 

The orbital-spin interplay due to the AL-VC is intu¬ 
itively understood in terms of the strong-coupling picture 
U IFband [37]: As shown in Fig. 2 (d), when the FO or¬ 
der rixz 2> riyz is realized, the nearest-neighbor exchange 
interaction has large anisotropy ^ Then, the 

stripe AFM order with Q = (tt, 0) appears if is not 
too small. Thus, the FO order/fluctuations and AFM 
order/fluctuations emerge cooperatively in the localized 
model, and such Kugel-Khomskii-type orbital-spin inter¬ 
play is explained by the AL-VC in the weak-coupling 
picture. 


LaFeAsO 



FIG. 2: (color online) (a) X^iQ) a-nd (b) Xx'^-y^iQ) obtained 
by the SC-VC method for LaFeAsO {z = 1). Note that 
other FO susceptibility (x'z^iQ)) and antiferro-orbital sus¬ 
ceptibilities iXxz/yziQ)) also develop secondarily, (c) Orbital 
(spin) enhancement factor Sc(s) as function of T for LaFeAsO 
(2 = 1) at r = 0.41. The charge and spin Weiss temperatures 
are 9c = 48 meV and 9s = 40 meV, respectively. Inset: 
Sc(s) as function of r at T = 50 meV. Note that U — 4.23r 
eV for LaFeAsO. (d) Orbital-spin interplay in the localized 
{dxz,dyz)-ovhit&\ model, known as the Kugel-Khomskii cou¬ 
pling. 

Next, we analyze the FeSe model, in which the ratio 
J/U is considerably small. In FeSe, the experimental 
mass-enhancement factor is ^ 10 for d^y-orbital, and 3 ^ 
4 for other d-orbitals according to the ARPES study [27]. 
Therefore, we put = z~^ for / 4 and z^^ = 3z~^ 

in the present study. We find that the peak of x®(g) 
moves from q = (tt, tt) to the experimental peak position 


q = (tt, 0) [4-7] for z^^ > 1.52“^, and the results of 
the SC-VC method are essentially unchanged for z^^ > 
1.52“^. Figures 3 (a) and (b) show the obtained x®(< 7 ) 
and x“ 2 _ 2 (g) for r = 0.25 [u = 1.76 eV) at T = 50 
meV in the case of 2 = 1. We see that only the FO 
susceptibility develop divergently [S'c « 50], whereas the 
AFM susceptibility remains small [Ss ~ 8], consistently 
with experiments for FeSe. The r-dependences of the 
Stoner enhancement factors at T = 50 meV are shown in 
the inset of Fig. 3 (c): With increasing r, Sc increases 
rapidly whereas Ss remains small. 

Figure 3 (c) shows the temperature dependences of the 
enhancement factors at r = 0.25. We stress that Sc 
approximately follow the Curie-Weiss behavior with the 
Weiss temperature 6c = 48 meV, which is consistent 
with the experimental Curie-Weiss behavior with positive 
6c in FeSe [2]. Since the spin Weiss temperature takes 
a large negative value {9s ^ — 20 meV), which is also 
consistent with experiments, one may consider that the 
orbital order in FeSe stems from causes other than spin 
fluctuations. 


IV. ORIGIN OF THE RELATION Sc > Ss IN 
FESE 

In this section, we discuss why the relation Sc ^ Ss 
{6c > 0 and 6s < 0) is realized in FeSe. First, we focus 
on the ratio between the Hund’s and Coulomb interac¬ 
tions J/U. It is intuitively obvious that the ratio J/U is 
an important control parameter for the orbital nematic- 
ity: For larger J/U, the local configuration of the two- 
electrons in the (dj,^, dy 2 )-orbitals is \dxz,\) <8) \dyz,'[), 
where the magnetic moment is Sj = 1 whereas the or¬ 
bital polarization is Uxz — nyz = 0. Thus, the smallness 
of J/U in FeSe is favorable for the emergence of the or¬ 
bital order without magnetization. 

Microscopically, as we discuss in Sec. II, the charge 
Stoner factor for Xx^-y2{0) is ac ~ {l-5J/U){l + UX‘^), 
where X‘^ is the charge AL-VC for orbital 2 or 3 at g = 0 . 
Since J/U = 0.0945 in FeSe, the orbital order is realized 
by relatively small AL-VC; A° ~ 0.9x°(0). In LaFeAsO, 
in contrast, large AL-VC of order ~ 2x°(0) is required to 
realize the orbital order. The obtained AL-VCs in both 
systems are shown in Fig. 9 (c) in Appendix D. 

We discuss why the AL-VC is important in the FeSe 
model with 6s < O'- As we explain in Appendix D. the 
T-dependence of the AL-VC is given as X'^ ^ X^TSs, 
where A is the three-point vertex that represents the in¬ 
terference between two short-living magnons. We find 
the relation A^ oc with a si I at low temperatures 
due to the good nesting between h-FSs and e-FSs [20]. 
Thanks to the strong enhancement of A at low temper¬ 
atures, the orbital order {ac = I) is realized even if 6s 
is negative. (Note that TSs decreases as T —0 when 
9s < 0.) Thus, serious diagrammatic analysis of the 
AL-VC is necessary to understand the rich normal-state 
phase diagrams in Fe-bases superconductors. 
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The enhancement of the nematic susceptibility due to 
the significant T-dependence of A^(oc T““) had been 
discussed in Refs. [20, 21, 51, 52]. However, the re¬ 
ported exponent a is not universal, since it depends on 
the bandstructure and temperature range. In Appendix 
E, we show the T-dependence of the three-point ver¬ 
tex for LaFeAsO and FeSe models for wide temperature 
range. It is found that a « 1 for T = 20z ^ I00z[meV], 
where z < I is the band-renormalization factor. Due to 
such large T-dependence of a, Xx^-y^iO) obtained by the 
present study follows the Curie-Weiss law only approxi¬ 
mately. 

Finally, we stress the importance of the orbital de¬ 
pendence of the spin fluctuation strength. Since the 
(ia;y-orbital h-FS (h-FS3) is absent in FeSe, the relation 
X2(3)(9) » xUq) ixtiq) = xIi;i,i{q)) is realized. This 
condition is favorable for the development of X° 2 _y 2 ( 0 ) 
since is enlarged by the spin fluctuations on the 

{dxz, dyz)-orbitals. More detailed analysis is given in Ap¬ 
pendix D. 


V. EFFECT OF THE MASS-ENHANCEMENT 
FACTOR 

Here, we study the effect of the mass-enhancement fac¬ 
tor: We study the FeSe model in the case of z~^ = 4 for 
dy 2 )-orbitals. The obtained Sc,s as functions of r 
are shown in the inset of Fig. 3 (d) at T = 12.5 meV. 
Here, Ss remains small even for r ^ 1 since the bare 
susceptibility is suppressed by z. In contrast. Sc is en¬ 
larged to 50 for r « 0.97, which is very close to the exact 
first-principles Hubbard model i7FeSe(?' = !)• The T- 
dependences of Sc,s are shown in Fig. 3 (d): Beautiful 
Curie-Weiss behavior with 9c = 12 meV is obtained by 
putting r = 0.97. 

To understand the similarity between the results in 
Fig. 3 (c) for z = 1 and the results in Fig. 3 (d) for 
z“^ = 4, we prove that both ac and as are independent 
of z under the rescaling T ^ zT and (t/, J) —)■ (t/, J)/z. 
Here, we assume that z^^ = z“^, and neglect the T- 
dependence of /i for simplicity. Under the scaling T —>■ 
zT, the Green function G(fc, n) at Matsubara integer n 
given in Eq. (4) is independent of z. For this reason, the 
bare susceptibility x°(9) = ^ G(fc-I-q, n)G(fc, n) 

is proportional to z. By following the same procedure, 
the three-point vertex A is scaled by z, and therefore the 
AL-VC X^{0) ~ TC/^ A(0; q)^X*(9)^ is proportional 
to z under the scaling T ^ zT and {U,J) —>■ {U,J)/z. 
Thus, both spin and charge irreducible susceptibilities are 
proportional to z, and both as and ac are unchanged 
under the rescaling T ^ zT and {U, J) {U, J)/z. That 
is, the Weiss temperatures 9s(c) are scaled by z. The 
validity of these scaling relations are confirmed by the 
numerical study in Fig. 3. 

It is possible to obtain z~^ by calculating the self¬ 
energy E(fc) together with self-consistently. In 

this case, fine tuning of r will be unnecessary since the 


FeSe 

Stoner enhancement factor 



Stoner enhancement factor 



r(meV) 


FIG. 3: (color online) (a) x*('7) and (b) obtained 

by the SC-VC method for FeSe (z = 1), for r = 0.25 at 
T = 50 meV. (c) T-dependences of the Stoner enhancement 
factors for FeSe (z = 1) at r = 0.25. Note that IJ = 7.21r 
for FeSe. Inset: Stoner enhancement factors as function of r 
for FeSe (z = 1) at T = 50 meV. (d) T-dependences of the 
enhancement factors for FeSe {z~^ = 4) at r = 0.97. The 
charge and spin Weiss temperatures are 9c ~ 12 meV and 
9s —7.5 meV, respectively. Inset: Stoner enhancement 
factors as function of r for FeSe (z“^ = 4) at T = 12.5 meV. 


relation as,c < 1 is assured if S(fc) and x*’‘^('7) are 
calculated self-consistently in two-dimensional systems 
(Mermin-Wagner theorem) [53]. This is our important 
future issue. 


VI. DISCUSSIONS 

A. Spin fluctuation strength and fe-dependent 
orbital-polarization below Tstr 

Here, we study the electronic states in the FO or¬ 
der Uxz 7 ^ Uyz established below the structure transi¬ 
tion temperature Tgtr, at which the shear modulus Gge 
reaches zero. According to the linear-response theory. 
Gee oc 1 - gXx2-y2(0), where Xx2-y2(0) oc (T - 9c)~^ 
is the electronic orbital susceptibility given by the SC- 
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VC theory, and g is the phonon-mediated Jahn-Teller 
energy [54], Therefore, Cee oc (T - Tstr)/{T - 9c), and 
Tstr = 9c + g is slightly higher than 9c due to the weak 
electron-phonon coupling (g « 10 50 K) [2, 17, 18]. 

Figure 4 (a) shows the T-dependence of Ss given by 
the RPA for LaFeAsO and FeSe for z = 1. Here, we 
introduce the orbital polarization —AEj2 (AA/2) for 
the c? 2 , 2 (j^ 2 )-level. We put Ss = 20 (5) for LaFeAsO 
(FeSe) at Tgtr = 50 meV, and assume a mean-field-type 
T-dependence; AE = ATotanh(1.74-\/Tstr/T — 1) with 
AEq = 80 meV. (For z~^ =4, the renormalized orbital 
polarization zAEq is just 20 meV.) In both LaFeAsO 
and FeSe, Ss are enhanced by AE, since as increases 
linearly with AE at q = (ti'jO) as discussed in Ref. [54]. 
In LaFeAsO, the magnetic order temperature T^ag in¬ 
creases from 9s to just below Tgtr since Ss is already 
large at Tgtr. In contrast, in FeSe, the enhancement of 
X®( 7 r, 0) is much moderate [55]. 

We also perform the self-consistent analysis of the 
orbital-polarization {AExz{k), AEyz{k)) and anisotropic 
X®(q), which is a natural extension of the SC-VC the¬ 
ory into the orbital ordered state [56]. The obtained Ss 
and fc-dependent orbital polarization are shown in Figs. 
4 (a) (inset) and (b), respectively. The parameters are 
r = 0.256 and 1 / 2:4 = 1.6. The difference An = n^z — riyz 
is 0.2%. The hole-pocket around F-point becomes ellip¬ 
soidal along the ky-axis due to the “sign-reversing orbital 
polarization”, in which AExz(0,k) — AEyz{k,0) shows 
the sign reversal as shown in Fig. 4 (c). Due to this sign 
reversal, Ss in the inset of Fig. 4 (a) tends to saturate 
below 40 meV [33]. Also, two Dirac-cone FSs appear 
around X-point when AEyz{TT,0) > 50 meV. These re¬ 
sults are essentially consistent with the recent ARPES 
studies reported in Refs. [27-33]. The obtained orbital- 
polarization {AExz{k),AEyz{k)) belongs to Big repre¬ 
sentation, and therefore it is consistent with the “d-wave 
orbital order” discovered in Ref. [31]. The d-wave orbital 
order is theoretically obtained by the mean-field approx¬ 
imation by introducing phenomenological long-range in¬ 
teraction [57], whose microscopic origin might be the AL- 
VC studied in this paper. 

In the present FeSe model with = 3, x'*(q, 0) has 
the maximum at q = (tt, 0 ), ( 0 , tt) without orbital order 
at T = 50meV, as shown in Fig. 3 (a). In the orbital 
ordered state, x'*(q, 0 ) at q = (tt, 0 ) increases as shown in 
Fig. 4 (d). These results are consistent with the neutron 
scattering study for FeSe for both T > Tstr and T < Tgtr 
[4-7]. Essentially similar results are obtained for zj"^ > 
1.5 at T = 50 meV. We verified that X®(q, 0) has clear 
maximum peak at q = (7r,0) even for zj”^ = 1.1 below 
T = 10 meV using 128 x 128 fe-mashes. Experimentally, 
Z 2 . 3 /Z 4 is about three [27], and the relation 22 , 3/24 ^1.3 
is predicted by the dynamical-mean-field-theory in Ref. 
[44]. 

According to Ref. [4], shows the broad max¬ 

imum at q = (tt, 0) at low-energies (w < 10 meV), and 
its strength is almost independent for T > Tstr- The 
magnitude of the low-energy spin susceptibility in FeSe 
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- 80 
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FIG. 4: (color online) (a) T-dependences of Ss for both 
LaFeAsO and FeSe models (z = 1). The FO order is in¬ 
troduced below Tstr = 50 meV. Inset: T-dependences of Ss 
for the FeSe model obtained by calculating the fe-dependent 
orbital polarization and X^(q) self-consistently. Ss tends to 
saturate below 40 meV due to the sign-reversing orbital polar¬ 
ization. (b) Self-consistent solution of the orbital polarization 
{AExz{k), AEyz{k)) in the orbital ordered state in the FeSe 
model at T = 50 meV. The shape of the C' 2 -symmetric FSs in 
(b) is consistent with the experimental reports [27-33]. We 
also show (c) the AE„.^(y,^i{k) along the fcj,( 3 ,)-axis, and (d) 
the C 2 -symmetric fo the orbital-ordered state. 


is one order of magnitude smaller than that in BaFe 2 As 2 
[58], whereas its magnitude would be comparable to that 
in LiFeAs [59]. This experimental report in FeSe will 
be consistent with the present theoretical result with the 
moderate Ss 10 in Figs. 3 (c) and (d). Note that ex¬ 
perimental dispersion relation in x*( 9 iw) for w < 100 
meV is qualitatively understood based on the present 
FeSe model by considering the band-renormalization fac¬ 
tor [7]. 


B. The ratio BsjQc for FeSe, NaFeAs, BaFe2As2 
and LaFeAsO as functions of JI'D 


In Fig. 5 (a), we summarize the ratio 9s 19c obtained 
in FeSe, NaFeAs, BaFe 2 As 2 , and LaFeAsO as function 
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FIG. 5: (color online) (a) Obtained QsjQc for four theoreti¬ 
cal models. Experimental values of Tmag/Tstr and ^NMu/Tstr 
are also shown. 0 nmr is the Weiss temperature of 1/TiT ob¬ 
tained by the Curie-Weiss fitting for T > Tgtr in LaFeAsO 
[60], BaFe 2 As 2 [61], NaFeAs [62]. In LaFeAsO, we derived 
0NMR « 95 K from the Curie-Weiss fitting of l/TiT above 

Tstr. In FeSe, we derived Tmag = —10 ~ —30 K from 

the Curie-Weiss fitting of 1/TiT below Tstr [2]. Thus, the¬ 
oretically expected relationships ^NMu/Tstr dsjOc and 

Tmag/Tstr > QsjQc are verified, (b) Obtained a’g for four 

theoretical models at ac = 0.98. We also plot for the 
hybrid models; LfilaFeAsO + (LaFeAsO TB model) and 
-ffpeSe + rRYii (FeSe TB model). 


of J/U. Numerical study for NaFeAs and BaFe 2 As 2 are 
presented in Appendix C. In NaFeAs and FeSe, in which 
J/U is smaller, the obtained Os!9c decreases to 0.4 and 
—0.4, respectively. In Fig. 5 (a), experimental values of 
T^mag/Fstr and ^NMu/Fstr are also shown, where 6»nmr 
is the Weiss temperature of 1/TiT above Tstr- Since 
Fstr = 6>c + ff (ff ~ 10 50 K) and 0 nmr = 9s, the rela¬ 

tion ^NMR/Tstr 9s 19c is expected theoretically. In ad¬ 
dition, the relation Tmag/Tstr ^ 9s/9c is expected since 
Tmag is substantially higher than 9s in the FO ordered 
state. These two theoretically predicted relations are ver¬ 
ified in Fig. 5 (a). Thus, the ratio 9s/9c is well scaled 
by the parameter J/U, consistently with the discussion 
in Sec. IV. 

Figure 5 (b) shows the critical value of the spin Stoner 
factor for ac ~ 1 in each model, a'g. It is found 
that Qfg' increases with J/U qualitatively. In addition, 
we plot a's foi' the FeSe (LaFeAsO) TB model with 
different Coulomb interactions: + rH^. 

In both FeSe and LaFeAsO TB models, monoton- 
ically increases with J/U, whereas is clearly small 


for FeSe TB model. There are two reasons why is 
smaller for the FeSe bandstructure. One reason is the 
absence of the d^y-orbital h-FS in FeSe: As we discussed 
in Sec. IV, the (ia;y-orbital spin fluctuations are un¬ 
necessary for the development of X^ 2 _y 2 ( 0 ) due to the 
AL-VC. Another reason is the smallness of the FSs in 
FeSe: We found numerically that a’g decreases when the 
size of the FSs is smaller, since the three-point vertex 
Am = 6xmiq)/^^Em, which is odd with respect to G, 
increases in magnitude when the particle-hole asymme¬ 
try is large: In fact, we analyzed the undoped LaFeAsO 
model with tiny FS pockets by introducing the posi¬ 
tive/negative potentials around the electron/hole FSs, 
and verihed that the orbital order is realized by small 
as- Recently, the advantage of the small FSs for the ne- 
maticity had been stressed by the renormalization group 
study in Ref. [63]. 


C. Summary 


The emergence of the electronic nematic order has 
attracted increasing attention as a fundamental phe¬ 
nomenon in strongly correlated metals. In this paper, 
we studied the origin of the nematicity in Fe-based su¬ 
perconductors, by paying the special attention to the 
nonmagnetic nematic order in FeSe. By applying the 
orbital-bspin fluctuation theory to the first-principles d- 
p Hubbard models, we succeeded in explaining the rich 
variety of the phase diagrams in Fe-based superconduc¬ 
tors, such as the nonmagnetic/magnetic nematic order in 
FeSe/LaFeAsO. The key model parameter to realize rich 
phase diagram is J/U; the ratio between the Hund’s and 
Coulomb interactions. In addition, the ratio 9s/9c tends 
to decrease as the size of the FSs shrinks, as discussed in 
Sec. VIB. 

In both FeSe and LaFeAsO, strong orbital suscepti¬ 
bility x° 2 _ 2 ( 0 ) oc (T — 9c)~^ with positive 9c is re¬ 
alized by tne strong orbital-spin interplay due to the 
strong-coupling effect, called the Aslamazov-Larkin ver¬ 
tex correction in the field theory. In the FeSe model, 
ferro-orbital order is established even when the spin 
Weiss temperature 9s is negative as shown in Fig. 3, 
since the three-point vertex (=the coupling between two- 
magnon and one-orbiton) increases at low temperatures 
as A ocT-O-5. In contrast, the spin-nematic susceptibility 
driven by the spin susceptibility should be T-independent 
if 0s < 0, as discussed in Ref. [7]. Therefore, we conclude 
that the nematicity in FeSe originates from the orbital 
order / fluctuations. 

The nematic orbital fluctuations might play important 
roles in the pairing mechanism in Fe-based superconduc¬ 
tors [64]. In FeSe, increases from 9 K to 40 K under 
pressure, accompanied by the enhancement of spin fluc¬ 
tuations [I]. At the same time, the system approaches 
to the orbital critical point since Tgtr decreases to zero 
under pressure. These facts indicate the important role 
of the spin-|-orbital fluctuations in FeSe. 
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Appendix A: Eight-orbital models for FeSe and 
LaFeAsO 


Here, we introduce the eight-orbital d-p models for 
FeSe and LaFeAsO analyzed in the main text. We first 
derived the first principles tight-binding models using the 
WIEN2k and WANNIER90 codes. Crystal structure pa¬ 
rameters of FeSe and LaFeAsO are given in Refs. [65] 
and [66], respectively. The obtained bandstructure and 
FSs in the LaFeAsO model are shown in Fig. 1 in the 
main text. In deriving the FeSe model, we introduce the 
fc-dependent shifts for orbital Z, 6Ei^ in order to obtain 
the experimentally observed FSs. In FeSe, we introduce 
the intra-orbital hopping parameters into JJpgSe in or¬ 
der to shift the dxy-orhiiaX band [(J 2 ,^/y^-orbital band] 
at (F, M, X) points by (0, -0.25, -bO.24) [(-0.24, 0, 
-1-0.12)] in unit eV. Such level shifts are introduced by 
the additional intra-orbital hopping integrals; Ji°n-site _ 
(5Er/4 -f 5Em/^ + 5Ex/2, = JAp/S - JAm/S, and 

j^nnn _ /16 J- SEm /16 — SEx /8. The bandstructure 

and FSs in the FeSe model are shown in Fig. 1 in the 
main text. 

We also explain the orbital-dependent Coulomb inter¬ 
action. The bare Coulomb interaction for the spin chan¬ 
nel in the main text is 


ini. 


I 2 5 ^ 3 ^' 


'Uh,h, 

J 2 5 

. 0 , 


h = h = h = h 
h = I3 ^ h = h 
/i = I2 7^ ^3 — ^4 
h = h h = I3 
otherwise. 


(Al) 


Also, the bare Coulomb interaction for the charge channel 
is 


-Uh,h, 

iniih.hh = + JhM 

J 2 1 

0 . 


h — h — h — h 

h = I 3 ^ h = li 
l\ — I2 ^ ^3 ~ ^4 

— ^4 7^ ^2 — ^3 
otherwise. 


(A2) 

Here, Ui^i, Uh, and Jip are the first principles Coulomb 
interaction terms given in Ref. [34]. 

Finally, we perform the band calculations for the or¬ 
thorhombic phase of FeSe and LaFeAsO, based on the 
experimental crystal structures. In both compounds, the 
obtained band splitting is too small to explain the large 
orbital polarization 60 meV) observed by ARPES 


studies. This result means that the orbital order orig¬ 
inate from the electron-electron correlation, which is not 
included in the band calculation. 

Figure 6 (a) is the non-magnetic bandstructure in 
the orthorhombic LaFeAsO obtained by the WIEN2k 
software. The spin-orbit interaction is not taken into 
account. The crystal structure parameters in the or¬ 
thorhombic phase is given in Ref. [66]. The orthorhombic 
structure deformation {a — h)/{a + b) is 0.3%. Due to the 
electron-phonon interaction, the four-fold symmetry of 
the bandstructure is slightly violated: The splitting be¬ 
tween the dxz- and dy^-bands, = Ey^ — Exz, is 

16 meV at X-point, and = 2 meV at F-point. 

Figure 6 (b) is the bandstructure in the orthorhom¬ 
bic FeSe. In the orthorhombic phase, the nearest Fe-Fe 
length is a = 2.6716A and b = 2.6610A, so {a — b)/{a + b) 
is 0.2% [65]. Here, the fe-dependent orbital shift to fit 
the ARPES bandstructure introduced above is not taken 
into account. In FeSe, = 14 meV at X-point, 

and = 3 meV at F-point. Thus, the sign re¬ 

versing orbital splitting observed in Ref. [33] cannot be 
explained by the band calculation. 

The splitting is reduced by the renormalization fac¬ 
tor z due to the self-energy. Since z ^ 1/3 in FeSe 
and LaFeAsO, the renormalized splitting at X-point is 
zAEband _ 5meV, which is one order of magnitude 
smaller than the experimental orbital splitting. There¬ 
fore, it is confirmed that the origin of the electronic ne¬ 
matic state in Fe-based superconductors is the electron- 
electron correlation. 


LaFeAsO 




FIG. 6: (color online) Bandstructure of (a) LaFeAsO and 
(b) FeSe in the experimental orthorhombic crystal structures 
obtained by the WIEN2k software. 
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Appendix B: Smallness of the VC for the spin 
susceptibility 

In the original SC-VC theory, the spin and charge sus¬ 
ceptibilities are calculated self-consistently, by including 
the MT-VC and AL-VC for the spin and charge suscep¬ 
tibilities [16, 49]. The strong orbital fluctuations are in¬ 
duced by the charge-channel AL-VC in Fe-based SCs, 
Ru-oxides and cuprate SCs [16, 37, 50]. In the main text, 
we studied the eight-orbital d-p Hubbard models based 
on the SC-VC theory, by taking the charge-channel AL- 
VC into account self-consistently. The obtained x®(q) 
is equivalent to the RPA since the spin-channel VCs are 
dropped. It is easy to verify that the charge- and spin- 
channel MT-VCs are negligible in the present model. 
However, the smallness of the spin-channel AL-VC is ver¬ 
ified only in the two-orbital Hubbard model in Ref. [50]. 

Here, we study the FeSe model using the SC-VC 
method, by taking the MT-VC and AL-VC for both spin- 
and charge-channels in order to confirm the validity of 
the numerical study in the main text. The charge (spin) 
susceptibilities are 

(HI) 

where + + The 

spin-channel AL-VC is given as 

p ar^h 

X [{yib-cAp + <i)yif-,g,h{-p) 

“b d)ye,f;g,h{ P')}-^m,m';c,d;g,hid^P) 

+ 2V'a,6;c,(i(P + 9)K 

where = ^c,h-m,g-d,m'{q\p) - 

^g,d-m,c-h,m'{.q\-p “ ?)■ Also, the expressions of 
the charge- and spin-channel MT-VCs are given in 
Ref. [49]. The double-counting second-order terms 
with respect to in -f i:AL.s(c) g^ould be 

subtracted [49] to obtain reliable results. 

Figures 7 (a) and (b) show the obtained X^io) and 
Xx^-y^iq), respectively, for U = 1.86 eV (r = 0.26) at 
T = 50 meV. The Stoner factors are obtained as as = 
0.907 and ac = 0.98. The obtained VCs 
and for c?a; 2 -orbital in the present self- 

consistent calculation are shown in Figs. 7 (c) and (d), 
respectively. Since is very small, the ob¬ 

tained charge and spin susceptibilities are very similar to 
the results in Fig.3 in the main text. Therefore, the valid¬ 
ity of the numerical results in the main text is confirmed 
by performing the very time-consuming self-consistent 
calculation with respect to and x®’'^(q). 



FIG. 7: (color online) (a) x'’(q) and (b) x^ 2 _j, 2 (q) given 
by the SC-VC theory, by calculating both spin- and charge- 
channel VCs self-consistently. The obtained results are quan¬ 
titatively equivalent to Fig. 3 in the main text. This fact 
means that the VC for the spin channel is negligible, (c) 
jj^MT+AL.s^^^ and (d) obtained by the present 

self-consistent calculation. 


Appendix C: Analysis of effective models of 
BaFe2As2 and NaFeAs 

In the main text, we introduced the first principles 
models for LaFeAsO and FeSe, and analyzed these mod¬ 
els by using the SC-VC method. Here, we also introduce 
the effective models for BaFe 2 As 2 and NaFeAs, and an¬ 
alyze them using the SC-VC method. 

In both BaFe 2 As 2 and NaFeAs, the FSs have rela¬ 
tively large three-dimensional characters. In addition, 
the unfolding of the bandstructure in BaFe 2 As 2 cannot 
be exactly performed because of its body-centered tetrag¬ 
onal crystal structure. Here, we introduce an simple 
effective BaFe 2 As 2 TB model L^BaFesAss magnifying 
the size of the d^jy-orbital hole-FS around k = (tt, it) in 
the LaFeAsO unfolded model, in order to reproduce the 
ARPES bandstructure in Bal22 compounds. Here, we 
shifted the da;y-orbital band at M point by -1-0.20 eV. As 
for NaFeAs, we just use LiLaFeAsO effective NaFeAs 
TB model, e.g., LI^aFeAs = ^LaFeAsO Considering that 
the FSs in NaFeAs in the = 0 plane are similar to 
the FSs in LaFeAsO. We use L^NaFeAs place of LI^FeAs 
given in Ref. [34]. 

The bandstructures and the FSs of the effective TB 
models for BaFe 2 As 2 and NaFeAs are shown in Figs. 
8 (a)-(c). Here, we perform the SC-VC analysis for 
the models Hm = Hm + (M=BaFe 2 As 2 , NaFeAs), 
where r(< 1) is the reduction parameter. We choose 
the parameter r to satisfy the charge Stoner factor is 
ac = 0.98; The obtained T-dependences of the spin and 
charge Stoner enhancement factors, Ss = (1 —as)“^ and 
Sc = (1 — ac)~^ respectively, are shown in Fig. 8 (d) 
and (e). As for BaFe 2 As 2 , both spin and orbital fluctu¬ 
ations strongly develop at T 50 meV in the case of 
r = 0.36. This result is consistent with experimental re- 
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Stoner enhancement factor Stoner enhancement factor 



T (meV) T (meV) 


FIG. 8: (color online) (a) Bandstructures of IFBaFe 2 Aa 2 
H^FeAs- (b) FSs of i/LFe 2 As 2 ^nd (c) FSs of//“.peAs- (d) T- 
dependences of the spin (charge) Stoner enhancement factors 
Ss(C) obtained in FfBaFejAsa- (Inset) The f7-dependences of 
the Stoner enhancement factors, (e) Spin and charge Stoner 
enhancement factors in //NaFeAs. 


lation Tinag ~ ’Tstr in BaFe 2 As 2 . As for NaFeAs, only 
orbital fluctuations strongly develop whereas spin fluc¬ 
tuations remain moderate at T ^ 50 meV in the case 
of r = 0.287. This result is consistent with experimen¬ 
tal results in NaFeAs [62], in which Tinag(= 40K) is more 
than ten Kelvin smaller than Tstr(= 53K). Thus, normal- 
state phase diagrams in BaFe 2 As 2 and NaFeAs are well 
explained by analyzing their effective Hamiltonians using 
the SC-VC method. 


Appendix D: Why are strong orbital fluctuations 
induced by tiny spin fluctuations in FeSe? 

In the main text, we studied the first-principles d-p 
Hubbard models for LaFeAsO and FeSe by applying the 
SC-VC theory. In both models, strong spin-fluctuation- 
driven orbital fluctuations are induced by AL-VC. In 
FeSe, we found that very small spin susceptibility x^ax 
is sufficient to realize the orbital order, consistently with 
experimental results. 

Here, we discuss why strong orbital fluctuations are 
induced by tiny spin fluctuations in FeSe. In Figs. 9 
(a) and (b), we show the spin and orbital susceptibilities, 
Xmax = X"(Q) and X“2_y2(0) =X2.2;2,2(g)+X3.3;3.3(9)- 
2 x 2 2-3 3 ( 0)1 ia the FeSe model and LaFeAsO model ob¬ 
tained by the SC-VC theory. Here, 32 x 32 fe-meshes and 
256 Matsubara frequencies are used. In both models, 
the charge Stoner factor is ac = 0.98 at T = 50 meV, 
and the obtained orbital susceptibilities show similar T- 
dependence. We set U = 1.76 {r = 0.25) in FeSe, and 
U = 1.74 (r = 0.41) in LaFeAsO, as we did in the main 
text. As for the spin susceptibility, in LaFeAsO, strong 
spin fluctuations develop at T = 50 meV {as = 0.98), 
consistently with previous theoretical studies [16, 37]. In 
FeSe, in contrast, Xmax is almost constant till T = 50 
meV {as = 0.87), consistently with experimental reports 
in FeSe. 

Now, we discuss why the spin fluctuation strength re¬ 
quired to realize ac ~ 1 is so different from LaFeAsO to 
FeSe. One reason is the difference in the ratio J/t/: Fig¬ 
ure 9 (c) shows the T-dependence of the AL-VC on d^z- 
orbital, A^'^’°(0) = A^ 2-’2 2 ( 0 ); obtained in the LaFeAsO 
and FeSe models. In both models, ac = 0.98 is satisfied 
at T = 50 meV. At T = 50 meV, the AL-VC for FeSe is 
about one-half of that in LaFeAsO. Thus, small AL-VC is 
enough to induce large orbital fluctuations in FeSe, since 
the charge Stoner factor is ac ~ (1 — 5J/U)U^2i^)- 

In Fig. 9 (d), we show that = 

^AL,c(o) _^AL,c. very Small for both FeSe and 

LaFeAsO. Here, “zero” represents the zero-Matsubara 
term (classical contribution) in Eq. (7) in Sec. H. Thus, 
non-zero Matsubara terms in the AL-VC are negligible 
in the present calculation (by chance). Note that the 
H^-term in AL-VC gives negative contribution. 

Another reason for the relation Xmax(I*'®^®) ^ 

Xmax(L3'FeAsO) at ac ~ 1 is the difference in the or¬ 
bital dependence of the spin fluctuation strength: The 
AL-VC for the sz-orbital is approximately given as 

X^^'%0) « 3rC/4^|A2,2;2.2;2.2(0;fc)|2y*(fc)2(Dl) 

k 

where we dropped the inter-orbital terms of x'* and A, 
and leave only the zero-Matsubara term in the Matsub¬ 
ara summation in Eq. (7) in Sec. H. If Eq. (Dl) is 
justified, only the spin fluctuations on dyz)-orbitals 
are important for the EO fluctuations. 

Eigure 10 (a) shows xKQ) EeSe and LaFeAsO for 
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FIG. 9: (color online) (a) Spin susceptibility x"{Q) (b) 

orbital susceptibilities for FeSe and LaFeAsO as 

functions of T. We put r = 0.25 for FeSe, and r = 0.41 for 
LaFeAsO. (c) and (d) for FeSe 

and LaFeAsO. 



FIG. 10: (color online) (a) xKQ), (b) C 2 = T,qXUQf, 
and (c) IA 2 P as functions of T in FeSe and LaFeAsO. 
(d) The approximate AL-VC for da; 2 -orbital = 

3i7'‘|A2(0; (0,7r))|2TC2. In both FeSe and LaFeAsO, the ob¬ 
tained X^^'^ qualitatively agrees to the exact numerical cal¬ 
culations in Fig. 9 (c). 


the same model parameters used in Fig. 9. As derived 
from Fig. 9 (a) and Fig. 10 (a), the ratio xV\Q)lx^{Q') 


is just 0.22 in LaFeAsO, whereas the ratio increases 
to 0.53 in FeSe, since the relation xliQ) ^ xKQ) 
ixliQ) X 2 iQ)) is satisfied in FeSe (LaFeAsO) because 
of the absence (presence) of h-FS3. This orbital depen¬ 
dence of the spin fluctuations in FeSe is favorable for 
realizing the FO fluctuations. 

To understand the mo del-dependence of the AL-VC in 
more detail, we calculate C 2 = show the 

result in Fig. 10 (b): The ratio cl^aFeAsO/j^FeSe jg j^g^ 
1.35 since the width of the peak of X 2 (^)^ around q ^ Q 
is much wider in FeSe. We also examine the square of the 
three-point vertex for (ia; 2 -orbital A 2 = ^ 2 , 2 ; 2 , 2 -, 2 , 2 {q, k) at 
q = 0 and A; = Q in Fig. 10 (c). In both models, the rela¬ 
tion IA 2 P oc with a « 1 is satished for wide tempera¬ 
ture range: Such strong T-dependence of the charge-spin 
coupling A 2 is essential for realizing the orbital fluctua¬ 
tions, so it should be taken into account in the numer¬ 
ical calculation. As results, we obtain a crude approxi¬ 
mation for the AL-VC, = 3C/^|A2(0; (0,7r))pTC'2, 

and show the result in Fig. 10 (d). This crude approxi¬ 
mation qualitatively reproduces the exact numerical re¬ 
sults for both FeSe and LaFeAsO given in Fig. 9 (c). 

In summary, in both LaFeAsO and FeSe, strong or¬ 
bital fluctuations are induced by AL-VC for the dxz{yz)- 
orbital, A^3j'^(0). In FeSe, very small spin susceptibility 
Xmax is sufficient to realize the spin-fluctuation-driven or¬ 
bital order, because of both the smallness of J/U and the 
largeness of € 2 - Strong T-dependence of A 2 is essential 
for realizing the orbital fluctuations due to AL-VC. 


Appendix E: Strong T-dependence of the 
three-point vertex 

In this paper, we found that the strong orbital fluc¬ 
tuations in Fe-based superconductors originate from the 
AL-VC for the orbital susceptibility. The moderate in¬ 
crement of the AL-VC at low temperatures shown in 
Fig. 9 (c) gives rive to the Curie-Weiss behavior of 
X°2_y2(0). For the increment of the AL-VC, the strong 
T-dependence of the three-point vertex, shown in Fig. 10 
(c), plays the significant role. Its strong T-dependence in 
Fe-based superconductors had been pointed out in Refs. 
[20, 21, 51, 52]. 

Here, we calculate the three-point vertex for LaFeAsO 
and FeSe models for wide temperature range with high 
numerical accuracy, using 512 x 512 fc-meshes and 2048 
Matsubara frequencies. Figure 11 shows the square 
of the three-point vertex for d^^z-orbital A 2 ( 0 ;( 5 ) = 
A2,2;2.2;2,2(0; Q) for T > 10 meV. In both LaFeAsO and 
FeSe models, the coefficient a of |A2(0;Q)p oc T“ de¬ 
pends on the temperature range. In both models, a « 1 
for T = 20 ^ lOOmeV, so the numerical result in Fig. 
10 (c) is confirmed by this accurate calculation. When 
the band renormalization due to z < 1 is considered, the 
relation a « 1 is realized for T = 20z ~ 1002;[meV]. 

For T < 20z [meV], |A2(0;(5)p saturates since the 
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FIG. 11: (color online) T-dependences of |A2(0;(5)|^ in 
LaFeAsO and FeSe models. 512 x 512 fc-meshes are used. 


temperature is smaller than the energy scale of the nest¬ 
ing. For T > lOOz [meV], the relation a « 2 is realized 
as discussed in Ref. [51, 52]. Note that the chemical 
potential p, becomes higher than the of the hole-band at 
r point when T is higher than lOOz (300z) [meV] in the 
FeSe (LaFeAsO) model. 

In Ref. [20], we reported the relation X^^’'^(0)/r oc 
ji-o.5(i_q,^)-i baged on approximate calculation for the 
five-orbital LaFeAsO TB model. The factor T~^-^ orig¬ 
inate from |A 2 |^. However, we performed more careful 
numerical analysis, and found that the approximate re¬ 
lation |A 2 |^ ~ T~^ (a ~ 1) is realized for T ~ 50^ [meV] 
in the five-orbital model. 


Appendix F: Two definitions of the averaged 
Coulomb and Hund’s interactions 

In the present study, the ratio between the intra-orbital 
Coulomb interaction and Hund’s interaction, J/U, is the 


essential control parameter. In this paper, we follow 
the Hubbard-Kanamori definition: U = ^ 

J = -^ J2i>m Ji,m- This definition is used in Ref. [34]. 
By using the Slater integrals [67], they are expressed as 
U = + ± (F2 + f 4), [/' = F’O - i (f2 + f 4), and 

J = ^ -I- F'^) [67]. According to the first-principles 

cRPA method [34], the ratio J/U is 0.0945 for FeSe, and 
the relation 17' = C7 — 2 J is approximately satisfied. 

In the first principles studies, another definition of the 
averaged interaction is used frequently: U = F^ and J = 
^ (F^ + F'^y This definition is used in Refs. [70, 71]. 

They are equivalent to U = (j2Li Ui,i + Y.i^m ^/.m) 

and J=Us-y-, {Uim - Ji,m) [68, 69]. 

Therefore, if we assume U' = U — 2J, which is actually 
satisfied well in Ref. [34], we obtain the relations U = 
U — ^J and J = I J. Thus, J/t7 = 0.224 obtained by the 
first principles study for FeSe in Ref. [70] corresponds to 
J/U = 0.127. Also, J/U = 0.294 for LaFeAsO obtained 
in Refs. [70, 71] corresponds to J/U = 0.157. 


As shown in Ref. 5 (b), the value of remains small 
(~ 0.9) in the FeSe TB model with [J/U = 

0.12) or with F^aFeAsO [j /U = 0.134). In each case, the 
obtained F-dependences of Ss and Sc are qualitatively 
similar to those shown in Fig. 3 (c). Therefore, the main 
results of the present study are unchanged even if J/U 
in FeSe is slightly larger than 0.1. 
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